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We discuss particle entanglement in systems of indistinguishable bosons and fermions, in finite 
Hilbert spaces, with focus on operational measures of quantum correlations. We show how to use 
von Neumann entropy, Negativity and entanglement witnesses in these cases, proving interesting 
relations. We obtain analytic expressions to quantify quantum correlations in homogeneous D- 
dimensional Hamiltonian models with certain symmetries. 
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The notion of entanglement, first noted by Einstein, Podolsky and Rosen [I], is considered one of the main features 
qq \ of quantum mechanics, and became a subject of great interest in the last years, due to its primordial role in Quantum 
Computation and Quantum Information [2h5|- Despite being widely studied in systems of distinguishable particles, 
less attention has been given to the case of indistinguishable ones. In this case the space of quantum states is restricted 
to symmetric (S) or antisymmetric (A) subspaces, depending on the bosonic or fermionic nature of the system. 

Entanglement of indistinguishable systems is much subtler than that of distinguishable ones, and there has been 
distinct approaches to its treatment, resulting in different notions like quantum correlations [6H8|; entanglement of 
■ modes and entanglement of particles [10]. In Zanardi's entanglement of modes @, as well as in Wiseman and 
Vaccaro's entanglement of particles, a Fock space is associated to the state space of a quantum system composed by 
several distinguishable modes, which allows one to employ all the tools commonly used in distinguishable quantum 
systems. In this work we will deal with the notion of quantum correlations which calls for different tools. 

Several notions of quantum correlations have been proposed in the literature, which agree in some respects, but 
differ in others. According to Eckert et al. Q, the pure states with no quantum correlations are those described 
. by a single Slater determinant for fermions, or a single Slater permanent formed out of a single one particle-state 
in the bosonic case. Li et al. [3 base their analysis on the resolution of the state in a direct-sum of single-particle 
states. Gihardi and Marinatto [tJ relate the notion of entanglement of quantum systems composed of two identical 
constituents to the impossibility of attributing a complete set of properties to both particles. It is important to 
note that these different definitions agree in the fermionic case, showing that the correlations generated by mere 
. anti-symmetrization of the state due to indistinguishability of their particles do not constitute truly as entanglement, 
or equivalently, states described by a single Slater determinant (Hartree-Fock approximation), which are eigenstates 
... of the free fermions Hamiltonian (single-particle Hamiltonian), have no quantum correlations. On the other hand, 
k** such definitions may disagree with each other in the bosonic case. Entanglement of indistinguishable fermions is far 
simpler than that of indistinguishable bosons. The definitions by Li et al. @ and Gihardi and Marinatto 0], although 
distinct, result in the same set of pure bosonic states without quantum correlations, which is greater than that defined 
by Eckert et al. @. Interestingly, as in the fermionic case, the former set corresponds to the eigenstates of the free 
bosons Hamiltonian, which is expected not to possess quantum correlations. 

Once one has opted for a certain notion of entanglement, the next step is to devise a method to calculate it. There are 
some interesting operational quantum correlation measures like the Slater concurrence @ for two fermions/bosons of 
dimension A('H 4 ®'H i )/S{'H 2 ®'rl' 2 )\ the von Neumann entropy of the single-particle reduced state for pure states of two 
particles [l,[ll|; the linear entropy of the single-particle reduced state of N-fermion pure states [l^- In a previous work 
[13| . we have shown how to calculate optimal entanglement witnesses for indistinguishable fermions, and introduced a 
new operational measure. With our witnesses we can calculate the generalized robustness of entanglement for systems 
with arbitrary number of fermions, with single-particle Hilbert space of arbitrary dimension. Interestingly, in the case 
of two fermions with a four-dimensional single-particle Hilbert space, the generalized robustness coincides with the 
Slater concurrence. All these measures have limitations, either conceptual or computational, and should be considered 
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complementary. The quantification of quantum correlations for general states, fermionic or bosonic, remains an open 
problem. 

In this work, as a natural extension of we will show how to calculate entanglement witnesses for the bosonic 
case, but they will not be optimal due to some subtleties of the uncorrelated bosonic states. We will show that 
functions of the purity of the single-particle reduced state quantify quantum correlations for pure states, with the 
caveat that for some special known values, the quantifier is inconclusive for bosons. This extends previous results by 
Paskauskas et. al and Plastino et. al [TH. We will also see that a simple shift in the well known negativity ( 
Neg(p) = \\p Ti Hj — constant) [l4| results in a quantifier of quantum correlations in bosons and fermions. Finally, in the 
context of entanglement in many-body systems [H, Il5l. Il6| . we will analyze homogeneous D-dimensional Hamiltonian 
models with certain symmetries, by means of the von Neumann entropy of the single-particle reduced state. 

This paper is organized as follows. In Sec|II]we consider quantum correlations in fermionic states, showing how the 
purity of the single-particle reduced state can be used as a measure for pure states, and the Negativity for the general 
case. In Sec lIIII the same analysis is made for bosons. In Sec II VI we discuss entanglement witnesses in bosonic systems. 
In Sec|V]we make some remarks about the different measures of quantum correlations, and discuss how they compare 
for pure states in the smallest dimension single-particle Hilbert space, proving some relations. In Sec I VII we show how 
to calculate quantum correlations in certain homogeneous D-dimensional Hamiltonian models. In the Appendix, we 
prove the expressions for the Negativity of bosons and fermions. We conclude in Sec lVIII 

II. FERMIONS 

Systems of indistinguishable particles have a more concise description in the second quantization formalism. The- 
refore we introduce operators with the following anti-commutation relations: 

{fl,f]} = {fijj} = 0, {fi,fj}=Sij. (1) 

fj and fi are the fermionic creation and annihilation operators, respectively, such that their application on the vacuum 
state (|0)) creates/annihilates a fermion in state "i". The vacuum state is defined such that that fi |0) = 0. 

As stated in the Introduction, the different definitions of quantum correlations agree with each other in the fermionic 
case, in the sense that the set of states without quantum correlations can be defined as follows. 

Fermionic state without quantum correlations: A fermionic state a 6 B(A(Hi<£) ■ ■ ■CS'H^)) has no quantum 
correlations if it can be decomposed as a convex combination of Slater determinants, namely, 

a = Y J Pic£ ■■■ 4|0><0|<--- 4, 5> = 1, (2) 

i i 

where a\. = Ylf—i u l kl fj ({a^ } is a set of orthonormal operators in the index k), U l is a unitary matrix of dimension 

dN , and {/, } is an orthonormal basis of fermionic creation operators for the space of a single fermion {Ti. d ). Note 
that uncorrelated pure states are single Slater determinants. The single-particle reduced states (oy (S!l) ) of a single 
Slater determinant have a particularly interesting form, and stand for the pure states in the "N-representable" reduced 
space (single-particle reduced space respective to the antisymmetric space of N fermions) jl~7j . 

Single-particle reduced fermionic state without quantum correlations: Given a pure fermionic state 
without quantum correlations, i.e. a single Slater determinant, \ip) = otat ...at^ |0), where {at.} are orthonor- 
mal, we have the equivalence: 

1 N 

a nsio = N^ a ^ l°>< l°** M =a \A,-A N 1°)' ( 3 ) 

i=l 

where oy (Sil) = Tri...Trjv-i(|'0) (V'l) is the single-particle reduced state (Tr-s is the partial trace over particle i). 
Therefore, if a is a mixed uncorrelated state, its single-particle reduced state in the "N-representable" reduced space 
is: 

oy = Tri...Tr N ^ t {a) = ^Pivl {Sll) - ( 4 ) 

i 

Now, aware of Eq.3, it is straightforward to conclude that shifted positive semidefinite functions of the purity of the 
single-particle reduced state can be used to measure the quantum correlation of a pure fermionic state, a result similar 



3 



to that obtained by Plastino et al. fjijj or Paskauskas et al. [llj. Using, for example, the von Neumann entropy 
(S(p) = Tr(-plogp)), we see that S(p r = Tr 1 ...Tr N _ 1 (\if;)(ip\)) > 5(er r(g( } ) = log AT, and thus a measure "E" for the 
quantum correlations of a pure fermionic state can be defined by a shifted von Neuman entropy of the single-particle 
reduced state. 

Shifted von Neuman entropy of entanglement for pure states: 

E(\TP)(i>\) = S(p r )-logN, (5) 

The case of pure states is easy due to the unique form of the uncorrelated single-particle reduced states (Eq.3), 
which is no longer the case for mixed states (Eq.4). Though not obvious, but straightforward to prove as we show in 
the Appendix, we can measure the quantum correlations of mixed fermionic states by the following shifted Negativity. 

Shifted Negativity: 

where Tj is the partial transpose over the i-th particle, and H-l^ is the trace-norm. If p is a single Slater determinant, 
its trace-norm is AT, and it is smaller in the case of a uncorrelated mixed state, as shown in the Appendix. Note, 
however, that we do not know if there are correlated fermionic states whose Negativity is null. 

III. BOSONS 

As in the previous section, we introduce operators to describe the bosonic system in the second quantization 
formalism. The operators satisfy the usual commutation relations: 
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Mi] = 0, b) =Sa, (7) 

where b\ and hi are the bosonic creation and annihilation operators, respectively, such that their application on the 
vacuum state (|0)) creates/annihilates a boson in state "i". The vacuum state is defined such that bi |0) = 0. 

As stated in the Introduction, the different notions of quantum correlations in bosons diverge from each other, 
resulting in two distinct sets of uncorrelated states. 

Bosonic pure state with no quantum correlations: A bosonic pure state \ip) <G S(Wf ® • • • ® %%), without 
quantum correlations, can be written as: 

N ° (bl 10) 

Definition 1. fofr> = J] n K ' (8) 
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Definition 2. \ip) = -^=(6^ |0) , (9) 



where 6^. = y ;fr _ T b\ is a set of orthonormal operators in the index i), U is a unitary matrix of dimension 

dN a , N Q is the number of distinct occupied states, and n c f >i is the number of bosons in the state fa. Uncorrelated 
mixed states are those that can be written as convex combinations of uncorrelated pure states. We clearly see that 
the set of states without quantum correlations according to Definition 1 includes the set derived from Definition 2, 
since the later is a particular case of the former, with N = 1. 

Definition 2 mirrors the case of distinguishable particles. Therefore one can use the entropy of the one-particle 
reduced state (S(p r )) and the usual Negativity (||/5 Ti || 1 — 1) to quantify the correlations. 

The problem is delicate for the Definition 1, for the equivalence between pure states without quantum correlations 
and the single-particle reduced states is no longer uniquely defined by the analogous of Eq.3. The shifted Negativity 
given by Eq.6 is still valid, but now we do know that there are correlated states with ||p :ri || 1 < AT. The entropy of 
the one-particle reduced state gives information about the quantum correlations, but as a quantifier it must be better 
understood. We know that a uncorrelated bosonic pure state, according to Eq.®, has the following one-particle 
reduced state: 

i=l 
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a r {4>i, <pj) is a matrix element of ay. The entropy of the one-particle reduced state assumes the special values: 

*K) = -£(^)i°g(^). (ii) 

1=1 

Note that < S(o~ r ) < log AT, and therefore when S(p r ) > log AT, the pure state p is quantum correlated. The pure 
state is also quantum correlated if S{p r ) is not one of the values given by Eq.ll. Take for example the case of two 
bosons: we have either N Q = 1, = 2 and thus S(a r ) = 0, or N = 2, = 1 and S(cr r ) = log 2. Given an arbitrary 
pure state p of two bosons, if S(p r ) = we can say with certainty that the state has no quantum correlations, but if 
S(p r ) = log 2 we cannot conclude anything, for either a state with no quantum correlations, e.g. \ip) = b^b^ |0), or 

a correlated one, e.g. \ip) = -^(ab^.b^ + Cjb^.b^. + Ckb^b^J |0), with a^ tk G and S(p r ) C (0,log3], could have 
the same von Neumann entropy for the one-particle reduced state. 



IV. WITNESSED ENTANGLEMENT 



In this section we present a bosonic entanglement witness. In one hand it is analogous to the fermionic entanglement 
witness we have introduced in a previous work [13j . but on the other hand it is not optimal, due to the complicated 
structure of the uncorrelated bosonic states. 

A Hermitian operator W is an entanglement witness for a given entangled quantum state p [Hjl, if its expectation 
value is negative for the particular entangled quantum state (Tr(Wp) < 0), while it is non-negative on the set of 
non-entangled states S (Ver e 6>, Tr(Wcr) > 0). We say that W op t is the optimal entanglement witnesses (OEW) for 
Pi if 

Tr(W optP ) = nun Tr(Wp), (12) 

where Ai represents a compact subset of the set of entanglement witnesses W. With OEWs we can quantify entan- 
glement (E(p)) by means of an appropriate choice of the set M. (l9| : 

E(p) =max(0,- min TriWp)). (13) 

In the fermionic case [l3[, restricting the witness operators to the antisymmetric space {W = AWA^}, the constraint 
{W < A} defines the Fermionic Generalized Robustness (R^)', while the constraint {Tr(W) = D a }, where D a is 
the antisymmetric N-particle Hilbert space dimension, defines the Fermionic Random Robustness (R^); and the 
constraint {Tr(W) < 1} defines the Fermionic Robustness of Entanglement ). These quantifiers correspond to 
the minimum value of s (s > 0), such that 

o = (14) 
1 + s K ' 

is a uncorrelated state (according to EqJ5]), where ip can be correlated or not in the case of R^, is uncorrelated in the 
case of R^, and is the maximally mixed state (A/D a ) in the case of Rf . 

The method for obtaining the OEW in the fermionic case is based on semidefinite programs (SDP) (20|, which can 
be solved efficiently with arbitrary accuracy. Now we will mimic the procedure for constructing W presented in [l^ . 
and try to obtain the Generalized Robustness for bosonic states. Consider the following SDP: 

minimize Tr{W p) 

d d d d 

£ - EE ••• E (clz\ ■••€sV.<-> 

ijv-l = l ii = l il=l jw-l = l 

subject tot , > 0, ( 15 \ 

Vcf G C, 1 < k < (N - 1), 1 < i < d, v ' 

SWS^ = W, 
W<S 7 

where d is the dimension of the single-particle Hilbert space, iS is the symmetrization operator, Wi N _^...i 1 j 1 ...j N _ 1 = 
bi N _ t • ■ ■ Wbj t ■ ■ ■ b^ N1 € B{Ji d ) is an operator acting on the space of one boson, and {bj} is an orthonormal basis 
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of bosonic creation operators. The notation W < S means that (S — W) > is a positive semidefinite operator. The 
optimal W obtained by this program is an entanglement witness, but it cannot be optimal, as we now discuss. 
For an arbitrary bosonic uncorrelated state a, the semi-positivity condition Tr(Wa) > is equivalent to: 

(0| b N b N _! • • • bi W b{ ■ ■ ■ 6] v „ 1 foj v |0) > 0, (16) 

for all orthonormal sets of creation operators {b\.}. This condition is taken into account in the second and third 
lines of Eq.15 by means of the semi-positivity of the operator &at_i ■ • -6i Wb\ ■ ■ ■ &jv-f Therefore, the entanglement 

witness W will not detect bosonic correlated states of the form b\ - ■ ■ b^ N _ 1 b\ r |0), where b N is not orthogonal to b\, a 
problem which does not arise in the fermionic case due to the Pauli exclusion principle. In numerical tests, we noticed 
that the quality of W improves with the increasing of the single-particle Hilbert space dimension. 



V. MEASURES INTERRELATIONS 



In this section we highlight the relationship among the measures of quantum correlations for fermionic and bosonic 
pure states in the smallest dimension, A(H 4 (&1-L 4 ) and S(H 2 (SiH 2 ), respectively. While the fermionic case resembles 
that of distinguishable qubits, the bosonic case is more intricate, due to the structure of the uncorrelated states. 

For pure states of distinguishable qubits, p = IV'XV'I G B(T-L 2 <g> H 2 ), it is well known the following equivalence for 
Generalized Robustness lZ g (p), Robustness of Entanglement 1Z e (p}, Random Robustness lZ r (p), Wooters Concurrence 
Cw(p), Negativity Neg(p), and Entropy of Entanglement E{p) |2ll - i24j : 

K g (p) = K e (p) = ±K r (p) = C w (p) = Neg(p) cx E(p). (17) 

Recall that E(p) is the Shannon entropy of the eigenvalues (A, 1 — A) of the reduced one-qubit state, and Cw = 

2y/\(l-\). 

For pure two-fermion states, p = \ip)(ip\ & B(A(H 4 ® H 4 )), we have found similar relations: 

nf{p) = Kf(p) = 2 -U^{p) - C£ (p) - \Neg{p) « E(p). (18) 

Note that Neg(p), and E(p) are the shifted measures. The relations between Robustness and Slater concurrence were 
observed numerically by means of optimal entanglement witnesses (l3| . and now we prove them. Based on the Slater 
decomposition — J^. ZiO ! 2i _ 1 a ! 2 i |0), where a\ = ^ fe Uikft, we can write the following optimal decomposition (viz 
Eq.14): 

Copt = Y~^_{p + t(j>opt), (19) 



'->opt 



-{a\4 |0)(0| a 3 ai + a\a\ |0)(0| a 4 a 2 ). (20) 



Now we show that when t = Cg(p), a op t is separable and in the border of the uncorrelated states. We know that 
the Slater concurrence of the state is invariant under unitary local symmetric maps <I>. We can always choose $ 
so that the single particle modes {aj} are mapped into the canonical modes {//} [25j]. Therefore <fr<r pt &' op t = 

i^dVWI + t&pt), where W) = J2* ZifLifL and <t>'o P t = WUl |0><0| hh + M |0)(0| / 4 / 2 ). 
The Slater concurrence of a' opt is given by Cg(cr' t ) = max(0, A4 — A3 — A2 — Ai), where {Aj}f =1 are the eigenvalues, 



in non-decreasing order, of the matrix J cr' opt a' opt , with a' opt = iJQAph)cr' opt {JQAph)\ being K, the complex conjugation 
operator, and U p h the particle-hole transformation. Consider the following matrix: 



yjv'optv'opt = y (1 + t) 2 (P'P' + Wo P t + VoptP') + t 2 <l>'opt<t>'opt)- (21) 

Note that u a' opt , p',(f)' opt " and their dual are all real matrices. With the aid of Eqs.19 and 20, it is easy to see that 
P'4>'opt = <P'o P tP' = 0, fioptfiopt = h<t>'o P u and that P'P' is orthogonal to 4>' opt $ opt . Thus Eq.21 reduces to: 



■opt" opt - / 1 + MVCf -r-^\/Yopt 



^p'p' + —.U'). (22) 
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The eigenvalues of y/ p'p' are easily obtained by means of its Slater decomposition, and the only non null eigenvalue is 
given by Cg(p'). ^J<j>' opt has just two non null eigenvalues, which are equal, given by (viz. Eq.20). Therefore the 
eigenvalues of the Eq|2T]are " {^s{p')i 1> IjO) "> an< ^ according to the definition of the Slater concurrence follows 
directly that (a' opt ) = if and only if t > Cg (p 1 ). 

We end this section by considering pure two-boson states, p = \ip}(i/)\ G B(S(H 2 <8> T~L 2 )) We have the following 
relations, which can be easily verified: 

Cl{p) = Neg(p) def . 2 (x E(p) def . 2 (23) 

In considering the measures corresponding to definition 1 of uncorrelated states (EqJ5]), we see that they are related 
differently, since the Negativity will always be zero for such states (Hp 7 *^ < 2). This is due the use of the upper 
limit in Eq |42l (viz Appendix). We could however, instead of using this upper limit, obtain analytically the values 
of /9' r< j corresponding to the uncorrelated pure states, which would be equal to p T< j = 1 or 2, and perform a 
similar analysis to that made for the S(p r )d e f.i in Eq llll Thus it would be possible to relate the Negativity and the 
Entropy of Entanglement according to definition 1. We see therefore that the relations between the distinct measures 
are similar to the distinguishable case when we consider the definition 2 (Eq|5]) of quantum correlations, possessing 
some discrepancies when we consider the definition 1. 



VI. HOMOGENEOUS D-DIMENSIONAL HAMILTONIAN 



In this section we see how to use the von Neumann entropy to quantify the quantum correlations in homogeneous 
D-dimensional Hamiltonian models, with the following properties: (1) the eigenstates are non-degenerate, and (2) 
the Hamiltonian commutes with the spin operator S z (thus S z and the Hamiltonian share the same eigenstates). 
Consider N particles of spin E, L D sites (with the closure boundary condition, L + 1 = 1), and an orthonormal 
basis {ct , C| } of creation and annihilation operators, where i = ..,«£>) is the spacial position vector, and (ct = 



— E, (— E + 1), (E — 1), E) is the spin in the direction S z 
(1) and (2), we have: 



If p = is one eigenstate according to the conditions 




(24) 
(25) 



where Eq. (|24[) follows from the translational invariance property of the quantum state due to the homogeneity of the 
Hamiltonian. while Eq.fJ^U) follows directly from condition (2). By condition (1) of non-degeneracy and the results of 
the previous sections, we known that the von Neumann entropy of the single-particle reduced state can be used as a 
quantifier of quantum correlations. Let us calculate it. 

We know that matrix elements of the reduced state are given by p r (ia,ja) = Tr(c\_c^ a \ip)(ip\) and, according to 

Eq. (|25|) . subspaces of the reduced state with different spin "ct" are disjoint. We can therefore diagonalize the reduced 
state in these subspaces separately. Eq. (|24|) together with the boundary condition fix the reduced state to a circulant 
matrix. More precisely for the unidimensional case (D = 1), given the subspace with spin "ct" and ■{c| cr }^_ 1 , the 
reduced state is given by the following L x L matrix: 



XL-l 



Pr = 



1 

N 



X2 



Xl 
XL-l 



■'■2 



XL-2 



Xl 
XQ 



XL-l 



XL-l\ 
XL-2 



Xl 
XQ J 



(26) 



XS = < C(k +S )aC' k!T >, 
Xq = 



L 



(27) 
(28) 
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where N a = X^j^i 71 ^- The terms x$ can be obtained by several methods, e.g. from two-point Green's function 
(one-particle Green's function). The eigenvalues {AJ}^ =1 of such circulant matrix are given by AJ = ^2^=o XkW j ' 
where Wj = cxp £^2, Thus the quantum correlations of that eigenstate can be calculated from S(p r ) = — AJ log AJ. 

For higher dimensions, given the subspace of a single-particle with spin "<r" and {ct }f =1 , the characteristic vector 
of its circulant matrix (e.g. the first line of matrix ) is given by, 

[D=2] : 

V c = ( [XOO ■ ■ -X(L-l)o] [%'"I(L-l)l] ••• 

••• [xo(l-i) •••a:(jr_i)(i_i)]) , (29) 

[D=3]: 

v c = ( «£> v™ ... v™^ ) , (30) 

where = ( [x 00 i ■ ■ ■ X(l-i)oz] [%mi • • • X{ L -i)ii\ • • • [x Q ( L -i)i ■ ■ ■ ^(l-i)(l-i)z] ) is the characteristic vector of the 
plane z = I, and xs x s y s z =< c \i mn ) cr c (i+s^(m+s y )(n+s,)a >• Thus, the eigenvalues {AJ}^ =1 of the reduced state are 
given by: 

L-l 

[D=2] : AJ = x im ^ +mL , (31) 

i,m=0 
L-l 

[D=3]:AJ= £ x imrl ^ +mL+ " L2 , (32) 

Z,m,n— 

where it>j = exp -rrr. If the eigenstate does not possess such properties, we can use the Negativity as a quantifier, 
but then we cannot provide analytic expressions. 

VII. CONCLUSION 

Entanglement between distinguishable particles is related to the notion of separability, i.e. the possibility of 
describing the system by a simple tensor product of individual states. In systems of indistinguishable particles, the 
symmetrization or antisymmetrization of the many-particle state eliminates the notion of separability, and the concept 
of entanglement of particles, referred in this work as quantum correlations, becomes subtler. If one is interested in the 
different modes (or configurations) the system of indistinguishable particles can assume, it is possible to use the same 
tools employed in systems of distinguishable particles to calculate the entanglement of modes. On the other hand, if 
one is interested in the genuine quantum correlations between particles, as discussed in the present work, one needs 
new tools. In this case, we have seen that quantum correlations in fermionic systems are simple, in the sense that the 
necessary tools are obtained by simply antisymmetrizing the distinguishable case, and one is led to the conclusion that 
uncorrelated fermionic systems are represented by convex combinations of Slater determinants. The bosonic case, 
however, does not follow straightforwardly by symmetrization of the distinguishable case. The possibility of multiple 
occupation implies that a many-particle state of Slater rank one in one basis can be of higher rank in another basis. 
This ambiguity reflects on the possibility of multiple values of the von Neumann entropy for the one-particle reduced 
state of a pure many-particle state. Aware of the subtleties of the bosonic case, we have proven that a shifted von 
Neumann entropy and a shifted Negativity can be used to quantify quantum correlations in systems of indistinguishable 
particles. Motivated by previous results with fermionic optimal entanglement witnesses, we have proven relations for 
robustness of entanglement and Slater concurrence for two-fermion systems with a four-dimensional single-particle 
Hilbert space, in particular showing that the Generalized Robustness and the Slater concurrence coincide for pure 
states. We have shown that the bosonic entanglement witness analogous to the fermionic entanglement witness is not 
optimal, due to the possibility of multiple-occupation in the former case. Nonetheless, numerical calculations have 
shown that the bosonic witness improves with the increase of the single-particle Hilbert space dimension. Finally, 
we have obtained analytic expressions for the calculation of quantum correlations in Homogeneous D-dimensional 
Hamiltonians. 
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Appendix: Negativity in fermionic/bosonic states 

In this appendix we calculate the trace-norm of the partial transpose of a uncorrelated fermionic/bosonic state, i.e. 

\\<J T '\\ = Tr[{<j T \a T * )3], thus proving the shifted negativity (Eq.6). We do so by the explicit diagonalization of the 

operator (a Ti a T i ). Consider first the case of a fermionic/bosonic pure state a = \ip) (ip\, as given by Eq.([2])/{8|), which 
can be rewritten as: 

a = C y^ e„e„rP w \(j>\<j>2—4>N) {4>N ■■■4>24>A Pit', (33) 

kit' 

with \tfj) = t-nP-n \(j>i4>2---4>N)i where (f>i, <pj are either equal or orthonormal, P^ are the permutation operators, 

e n is the permutation parity (e = ±1 for fermions, e = 1 for bosons), and C = (A!) -1 for fermions or C — 
[TV! OiI°i( n </>i')] _1 f° r bosons. From now on we omit the normalization C and introduce the following notation: 

P K \<h....<t> N ) = \-k{4>i..4n)) = k{<kM<h)-MM) . (34) 
Now we make the partial transpose on the first particle explicit: 

cr Tl = ^e^ev K'(^i)i"(02— ^jv)} (tt'C^at— </!>2)7r(0i)| ; (35) 



(a T ^=a T >; (36) 



^ e^'C^Cff/ W((f>i)^(h— 4>n)) 

7r,7r' ,7r,7r' 

(7r'(^jv...</>2)7r(0i)|7r'(0i)7r(02...0iv)) {n 1 {4> N ...4> 2 Wi)\ ; (37) 

7r' ,7T 

5^ W (7r(^i)|vf'(^i)) |tt(0 2 ...^)) (^(^...0 2 )| . (38) 

7T,7r' 

We analyze only the bosonic case, and the fermions follow by setting N = N and = 1. Consider the first 
line of Eql38l As states 4>i are not necessarily orthogonal, and may be the same, we have contributions when 
the permutations 7r',7r are equal and in some cases even when they are different. It can be seen that there are 
rik[{N — 1)!] permutations such that 7r'(</>i) = and for each of these there are Oi^i permutations tt such that 
tt(4>i) — <f>ki resulting in non null contributions (tt' '(4>n ■■■4'2)\ : k(4 i 2---4 i n)) ^ 0. If 7r(^>i) 4>k then the contribution 
is null (7r'(07v...(/)2)|7r(02---0A r )) = (simply note that the set {7r(</>2 ■•■</' jv)} always has states "</>*/', whereas 
{ir'{4>N---4'2)} has only n& — 1). The first line of Eq I3"51 thus reduces to: 

5> fc [(tf-l)!][I](n*!)] (39) 

k=l i=l 



Now we analyze the second line of Eq|3S] This term has non null contributions only if n((j>i) = n'((j>i). For permutati- 

(Af-l)! _ n^JN-l)'. 

(n k -l)\ n (»«,!) flWiO 



ons of the type tt(4>i) = tt'(4>i) = 4>k, the matrix \n(<j>2-.-4>N)) (k'(4>n---4>2)\ can assume ^ N J~^' = — 
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distinct combinations from the elements of the set {7r(02---0w)}- Note that there are permutations of type 

7r((/>i) = (j>k generating the same "ket" |7r(^2---^Jv)) ( or "bra" (7f'(<^iv-..^2)|)- Thus we have, 



n 



(40) 



%=i 



where \ipk) = Ej | 7r fc(02---0Ar)), being ir l k {(f>2---4>N) all the possible permutations such that it' l k ((j)i) = (j)k, and 



7T k {(f>2—<PN 

reduced to: 



7r£ (</>2 ■•■</' jv) ) = 5ij. We have then (ipkli'k') = 7^7— rr^fcfc' , and finally the second line of Eg 1551 is 



N N 



No 



N„ 



(41) 



k—1 i—1 i—1 k—1 

From Eq. (|39p and Eq.(|41l ) and remembering to reintroduce the normalization constant C, we obtain: 



N 



T A 



(E V^hI 



k=l 



N 



< N. 



(42) 



N No 

The last step follows by noting that ^ ~ N, and thus \/^k < N. As the trace-norm is a convex function, we 

fe=i fe=i 



can write for uncorrelated mixed states: 



E ; /',< , <E,Pj 



and we are done. 
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